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Abstract 


A carbon-nitrogen model has been constructed for some soil processes 
taking place in the FH-layer of a pine forest. It includes processes 
as decomposition, respiration, mineralisation, immobilisation and 
decomposers biomass formation. Most of the data used to construct the 
model was taken from the literature and a few from own measurements. 
A number of runs have been carried out. One of them is defined as the 
standard run. The other deal with the sensitivity of the model to 
changes in chosen parameter and to changes in structure. Mineralisation, 
in particular, shows to be very sensitive to changes in parameters 
describing the kinetics of decomposers' growth and to changes in the 
model structure. 


1. Introduction 


The study of soil processes in the project has been formalised in a soil 
nutrient model (named the N-model - present generation designated N 1). 

The overall aim of the model is to describe the competition for nitrogen 
(other nutriénts are of secondary priority) between roots, micrcorganisms, 
and soil (absorption) in an undisturbed forest soil ecosystem. For practical 
reasons, the model.has been divided into an organic (NIORG) and an inorganic 
(NIING) submodel. This development paper deals with the organic submodel 
(N10ORG). à; 


Several attempts have been made at constructing simulation models of the de- 
composition process. Many difficulties arise when, for example, the decomposi- 
tion processes, biomass formation and nitrogen mineralisation are to be organ- 
ized in a simulation. model. One requirement is that the model should have bio- 
logical significance, and another is that data obtained should fit into it. 
Lack of suitable data often limits the possibilities of constructing and vali- 
dating a model. For example, the processes of decomposition, mineralisation, 
biomass formation as well as the respective compartments often have to be 
studied on a small test site and a high degree of orgànisation in the studies 
is needed. 


Many investigations have been carried out on decomposition determined as 
weight loss. Less work has been done taking into account the actual changes 
taking place in the substrate during decomposition. It is well known, though, 
that the C/N quotient in substrates with little nitrogen decreases during the 
decomposition process (Smith & Douglas, 1971; Berg, to be published). The 
inorganic nitrogen formed is believed to be immobilised during the microbial 
breakdown of the substrate. The nitrogen is thus largely recycled in the 
substrate while carbon is released in the form of C05. The measured increas- 
ing percentage of nitrogen could be due to the presence of microorganisms 

and the substrate itself could have a lower rate of increase in the C/N 
quotient (Harmsen & Kohlenbrander, 1965). 


The present paper is intended to be a presentation of the state of the work 
and may be regarded as a working document for the present refinements of 
the model. It is intended that this paper should be followed by more com- 
plete descriptions as improvements and refinements are worked out. 


2. Objectives 


The objective of the organic part of the soil nutrient model (NIORG)are to 
describe the mineralisation and immobilisation of nitrogen with the 

time resolution of one day. The nitrogen mineralisation is an input to 

the inorganic part of the model (NTING). 


The model presented is modified from a submodel of decomposition devel- 
“oped in the US Desert Biome (Parnas & Redford, 1974). Although orig- 
inally developed for a desert, its general applicability is such that it 
was possible to validate it with data from the Norwegian IBP Tundra 
Biome studies. These data are quite different from those by which the 
model was constructed, but nevertheless gave acceptable results 
(Wielgolaski, Redford & Valentine, 1974). 


Our first aim is to make the presented model work properly with avail- 
able data from the FH-layer in a 120-year-old pine stand (Ih VA at 
Ivantjárnsheden). The soil is divided into three layers viz., the S- 
layer, a combined FH-layer, and the mineral soil down to 20 cm depth. 
The three layers behave differently in many respects. We preferred, how- 
ever to work with one soil layer in a first attempt to make the model 
work. A number of runs were carried out with the four available data of 
our own. It is our intention to run the three soil layers in a forth- 
coming N 1 model. This will be built around a number of modules. The 
present NIORG model constitutes the basis for the organic module of N 1. 


3. Verbal and graphical description of the model 


For each soil layer, a number of different structures can be chosen to 
build the model NIORG (see section 5). In our runs with the FH-layer we 
have been particulary interested in two of these structures and they 

are shown in figures 1-а and 1-b. With the structure shown in figure 1-а 
three separate decomposing substrates (organic materials is used sometimes 
in this paper as synonymous) are present in the layer, i.e. litter (L), 
dead organisms (DO) and humus (H). The second structure with only two 
substrates is shown in figure l-b. Here, humus and dead organisms have 
been combined in one decomposing substrate that we have called, 
provisoriously, soil organic matter (SOM or S) (see 4.5). Regardless of 


the structure, for each substrate three compartments (state variables) 


are defined (see 4.2): Ү(1) is protein carbon, Y(2) is other carbon and 
Y(3) is protein nitrogen in the substrate. Compartments X02,(decomposers 
biomass) and X01, (inorganic nitrogen in the soil) are common to the 

soil layer and independent of substrate concentrations. Inorganic nitrogen, 
(ХОЛ), can be present partially in soil water phase and partially adsorbed 

to soil particles. These amounts are represented by functions Gl and G52 
respectively (see 4.12 and 5.1.1). 


At each time interval (1 day in the model) the processes shown in 
figures l-a and l-b will modify the values taken by the compartments at 
the beginning of such time interval. An enlargement of the processes 
accompanying the decomposition of a given substrate is given in figure 1-с. 
This figure shows in a clear way the processes associated to carbon and 
nitrogen dynamics. The rate of processes are represented by the G-functions 
showed and explained in the same figure. The nomenclature used in 

figure 1-с for the G-functions and for the substrate compartments is of 
a general nature and refers to any substrate at any soil layer. Further 
on, in order to refer to a given G-functiom or compartment for one 
particular substrate we shall add to the G @r compartment symbol the 
capital letter representing the substrate. This means that GI9L, 61900 
and GI9H, for example, will represent mineralisation from litter, dead 
organisms and humus respectively. In some cases and in order to simplify 
writing, the vectorial notations has also been used; in such a case, 
Gl9(m) will mean the mineralisation from substrate m. Once the rates 

for each substrate have been calculated the total rate for the soil 
layer can be calculated by adding the rates for each substrate. So, 

the total nitrogen mineralisation for the ЕН! layer is represented by 
GI9T = GISL + 61900 + 619Н, or, in vectorial motation by: GI9T = eal im). 
The method of calculating the rate of processes for any substrate 

will be discussed below. However, two variables of interest must 

first be presented: 


- the carbon-nitrogen ratio of the substratie, defined at each time t, 
as the sum of the carbon compartments, Y (1) and Y (2), divided by 
the organic nitrogen compartment Y (3) (see 4.2 and 5.1.4) and 


- the critical carbon-nitrogen ratio, which does not depend on time t, 
and is defined as: average carbon fraction in decomposer cell div- 
ided by the product of average nitrogen fraction in decomposer cell 
and . decomposer assimilation efficiency ((see 5.1.6). 


Since practically all processes shown in figure 1-с depend on G 13 
(biomass formation) at each time t this process is calculated first. Some 
of the factors on which G 13 depends are:Y (1), Y (2), Y (3), NH; > decom- 
poser biomass, temperature, moisture and the optimal growth for the sub- 
strate. Thus, in order to calculate G 13, the temperature and moisture . 
conditions at time t must first be defined. This is done by means of the 
driving functions (see 5.2). G 13 can be calculated as the "stationary 
growth" (Parnas, 1975a) or as the "mean growth" (Parnas, 1975b). In any 
run of the model, one of both alternatives must be chosen, and this is 

done by attaching a convenient value to one of the parameters of 

the model (see 5.1.10). Once G 13 has been calculated, the decomposer 
demand on carbon and nitrogen can be calculated as P018, G13 and 

Р020. 613 respectively, where P018 and P020 are two parameters giving 

the carbon and nitrogen concentration in the decomposer cell (see fig. 1-с). 
The. total carbon decomposition, G 14, is calculated by dividing this 
carbon demand by the assimilation efficiency (5.1.12). On the other hand, 
G 14 minus the carbon demand gives the respiration rate,G 17 (5.1.16). 


Nitrogen fixation, G 23 is calculated by multiplying the nitrogen de- 
mand by the fraction of decomposer biomass which fixes nitrogen (5.1.13). 
The amount provided by G 23 is then subtracted from the total nitrogen 
demand; the remaining amount is multiplied by the inorganic nitrogen 
fraction, NHj/(Y (3) + ИНД), and in this way immobilisation, б 20, 

is calculated (5.1.19). The sum of G 20 and G 23 is subtracted from the 
total nitrogen demand and in this way we obtain the amount of nitrogen 
demand provided by Y (3), i.e., the organic nitrogen demand (arrow with- 
out number in fig. 1-с). 


We then calculate the protein-carbon decomposition. This is done as follows 
(5.1. 149: 


- if the carbon-nitrogen ratic of the substrate is greater 
than the critical carbon-nitrogen ratio, G 15 equals the organic nitro- 
gen demand times the carbon-nitrogen ratio of the protein. 


- if the opposite happens, G 15 is calculated as G 14 times the protein- 
-carbon fraction, Y (1)/(Y (1) + Y (2)). 


Other carbon decomposition is calculated as G 14 minus G 15 (5.1.15). Pro- 


tein-nitrogen decomposition, G 18, is calculated as G 15 divided by the 


carbon-nitrogen ratio of the protein (5.1.17). Mineralisation, G 19, is 
calculated by subtracting the organic nitrogen demand from G 18. It can 
be seen that if protein-carbon has decomposed according to the first way 
defined above, G 19 will equal zero. 


The remaining processes shown in fig. 1-с are formulated in a very simpli- 
fied way by the model (see 4.11 and 5.1.11 for mortality, 612; 4.7.3 and 
5.1.23 for external breakdown, G39, and 5.1.25 for nitrogen uptake by 
roots, G70). Regarding the inputs also shown in the figure, they 

clearly depend on the substrate in consideration (see figs 1-а and 
l-b). If the substrate is litter, the inputs will represent the different 
chemical fractions in litter formation, a driving function (see 5.2). 

If the substrateis DO (dead organisms), the inputs are given by 612, 
mortality, times the fraction of the given chemical compound in the 
decomposer biomass, plus the amount coming from other soil layers 

(see 4.6 and 5.1.24). In case the substrate is humus, the inputs will 
represent external breakdown coming from litter and dead organisms. 
Finally, if the substrate is SOM (soil organic matter), inputs will 
represent mortality plus transport plus breakdown (see 4,5 and 5). 


4. General assumptions and definitions 


In order to formulate the mathematical descriptions for the processes it 
was necessary to make a number of assumptions for formulating equations of 
the processes. One assumption is that the litter or organic matter being 
decomposed is homogeneous. More specific are the following assumptions/ 


simplifications taken from Parnas (1975a). 
Carbon loss from the organic substance under decomposition 15 pro- 


portional to the growth rate and requirements for maintenance en- 
ergy for the microorganisms decomposing it. 


- The growth rate of the microorganisms under certain specific en- 
vironmental conditions (temperature, moisture, pH,etc.) and a 
given type of litter or other organic substance increase 
with increasing carbon and nitrogen concentrations. 


- The decomposing microorganisms use the nitrogen compounds as 
nitrogen sources according to their relative availability. 


- The release of nitrogen from carbon-nitrogen compounds (C-N com- 
pounds) in the organic substance under decomposition is pro- 
portional to the loss of carbon from the same C-N compounds. 
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- As energy and carbon source, the decomposing microorganisms can use 
carbon bound or not bound to nitrogen. 


4.1 The soil layer (FH-layer) 

The FH-layer extends from the S-layer down to the middle of the 
bleach soil. The only litter appearing in this layer is root litter; 
and the amount is given in Table 1 (Persson, pers. comm.). There is 
an input of solid material from the S-layer consisting of humus and 
dead organisms. The thickness of the layer varies at the site stud- 
ied, but this variation is not considered here. The amount of soil 
organic matter was estimated at 1300 9 m? (Staaf & Berg, 

1976). Charcoal from a forest fire was found at the site but 

is neglected here. 


4.2 Substrate quality in the model and C/N quotients 

The substrates consist of various chemical components in differ- 

ing amounts. It is believed that these have different rates of de- 
gradation and it seems that lignin is degraded more slowly than cellu- 
lose and hemicelluloses,which are degraded more slowly than water sol- 
uble substances (Berg, to be published). In the NiCRG model the 
substrate quality is so far used in the following way: The model 
subdivides each of the substrates viz., 


- Litter, dead organisms and humus into three subgroups on a chemical 
basis. The easily available carbon viz. that in proteins and nu- 
cleic acids, makes up the Y(1)-compartment. It is defined as "the 
amount of carbon bound to nitrogen". 


- Carbon compounds not containing nitrogen, like cellulose and hemi- 
cellulose and - so far - C in lignin (in spite of "lignin 
nitrogen"), constitute the subcompartment Y (2). 


- The third compartment Y (3) contains the amount of nitrogen bound 
to carbon Viz., the amount bound to proteins and nucleic acids. 
These C and N compounds are used to define a C/N ratio only taking 
into account the protein-nitrogen of the substrate as nitrogen 
source. The quotient is referred to as C/N). Since the N4-nitrogen 


n 


cannot be measured in substrates where microorganisms are present, it 
must be calculated. In the model outputs we can obtain control values 
where the actually measured C/N is given. This quotient is referred 
to as measured C/N and corresponds to the traditionally used C/N 
quotient which includes microorganisms. 


The reason for this subdivision and the special C/N4-quotient is the 
different conditionsfor mineralisation of nitrogen and the utilis- 
ation of carbon. It is uncertain whether this subdivision is the best 
one for our system, and work is being carried out onthe possibility of 
including another compartment for nitrogen viz.lignin-bound nitrogen. 
There is now no difference between e.g. lignin and hemicelluloses in 
decomposition rate in the model and all substances classified as 
Y2)-compounds in litter have the same decomposition rate. 


4.3 Litter 
Litter is found in the three soil layers but its components vary ac- 
cording to the layer. 


The FH-layer litter is made up of root litter exclusively viz., 
heather roots, pine roots and lingonberry roots. These litter 
components, as well as those in the S-layer, have different decompo- 
sition rates. The different rates depend on the type of litter апа 
the diameter of the roots, e.g. pine root litter is decomposed twice 
as fast as heather roots, and litter of 3 mm diameter about twice as 
fast as that of 6 mm (Berg, to be published). An average decomposi- 
tion rate is used in the model. 


4.3.1 Litter inputs 

There are very specific abiotic factors, such as temperature and moisture, 
that govern the formation of litter. However, owing to lack of information, 
it has so far been necessary in the model to give the litter inputs at 
selected times of the year independent of abiotic factors. 


The litter input to the FH-layer consists only of root litter, viz., pine 
roots of diameters less than 5 mm, heather and lingonberry roots. The 
present assumption is that two litter inputs equal in size take place 
annually, during November and April. 
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4.4 Input of solid organic material from the S-layer 

The litter that falls on the S-layer or is formed there is decomposed to 
such an extent that it is transformed into soil organic matter (SOM) or 
humus already in the S-layer (Staaf & Berg, in manuscript). This means 
that the lowest part of the S-layer is made up of SOM and the input to the 
FH-layer thus becomes either SOM or humus and dead organisms (see the 
discussion in section 4.5 "SOM versus humus and dead organisms"). 


The mechanisms for downward transport of litter through the moss/lichen 
layer are not known. The transport may, for example, be influenced by the 
snow, pressing down the litter particles in the upper parts of the moss 
layer. Moss growth may also have an effect. As the mosses grow, the needles 
descend in the profile in their original position. For the purpose of the 
model we have made the following approximation of the transport: the SOM 
(humus and dead organisms) are transported at a constant rate during 

the time when temperature of the layer is above 00C. 


The transport time through the S-layer determines the amount of carbon that 
enters the FH-layer owing to the release of C02. The nitrogen is not de- 
pendent on transport time since it appears that none or very little nitro- 
gen is given off from the litter particles during breakdown (Berg, to be 
published). 


Measurements and calculations of times for the downward transport of needle 
litter through the S-layer gave a rough figure of 5-6 years (Berg 1974; 
Berg, in manuscript) at which the litter might be decomposed by about 90%. 
Using a value of 200 g · m^? yr7l litter input to the S-layer (Flower- 
Ellis, pers. comm.) we would thus obtain 20 g * m"? - yr-l of input to 

the FH-layer. This quantity of 20 g thus retains all the nitrogen present 
in the original 200 g of litter (0.42). 


4.5 Soil organic matter (SOM) versus humus and dead organisms 

Litter under decomposition finally reaches a stage where it can no longer 

be described as litter. At this stage the word humus is applicable only as 
unidentifiable substances remain. The term soil organie matter (SOM) is as 
applicable and has been used much in modelling work to designate the 
Structurally amorphous mass of humus, among other things, without specify- 
ing what is included. The SOM includes humic substances, litter remains, and 


dead microorganisms. Considering the average amounts of living microorganisms 
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апа their high nitrogen content, their weight and amount of nitrogenare 
sometimes subtracted when calculating the amount of SOM. There are different 
ways of regarding the visually unidentifiable organic substance present in 
the FH-layer. It can be regarded as a rather homogeneous mass or as con- 
sisting of different components. In the present work and in the N10RG model- 
ling we chose to distinguish carefully between what we call soil organ- 

ic matter and what we call humus: Among the reasons for this are diff- 

erent approaches to the nitrogen mineralisation. The meanings of the 
different terminologies used here are discussed below. The model was 

run with the two following approaches: 


- As soil organic matter 
Soil organic matter is taken to mean the amorphous substance including 
the nonidentifiable humus substance and dead organisms. The figure for 
the living organisms - in reality the fumgal and bacterial biomasses - 
has been subtracted. In this case the SOM contains protein-nitrogen 
and protein-carbon as well as "other саг%оп". The SOM has a C/N quot- 
ient of about 50 (Staaf & Berg, in manuscript). 


- As humus and dead organisms 
Іп this case the amount of dead organisms in the SOM is calculated and 
subtracted. We thus obtaintwovery different substances. The humus con- 
tains no easily degradable carbon-nitrogen compounds and the protein- 
nitrogen subcompartment is therefore set at a low value and so is the amount of 
protein-carbon. The "other" carbon subcompartment has a rather high 
value, and the humus compartment thus has a very high value for the 
С / Ny quotient. The slowly available mitrogen in the humus has not 
been taken into account so far. 


In the dead organisms compartment the C / № quotient is 11.3 
(Söderström, pers. comm. ). 


4.6 Nitrogen inputs 
The main nitrogen input appears to come witth the litter formation 
(see 4.3.1) and might be 0.46 g- m ^ - yri. 


The other large nitrogen input probably takes place vie the transport of SOM 
(or dead organisms) and is thus protein nitmnogen from the е adir 
input,which might reach an order of magnitude of 0.1 д · m^ yrs iu 


discussed in section 4.4. 
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The nitrogen fixation in the FH-layer is very small, probably of an ог- 
der of magnitude of 0.01 g · me - yr! (Granhall, pers. comm.). It has 
been included in the model, however, and it is assumed that the nitrogen 
fixation is carried out by 2% of the biomass. The nitrogen fixation has 
the same temperature and moisture dependence as growth and decomposition 
processes. 


Dry and wet deposited N (мна) goes to the S-layer and is assumed to be 
further transported by water or taken up by the decomposer organisms and 
the vegetation. It has thus been incorporated in other processes, and is 
not a separate process in the present model. The faeces fall is rather 
small but might be included in the model later on. No denitrification 
appears to be taking place at the site and the process has thus 

not been included in the model. 


4.7 Decomposition 


General 

The decomposition process depends, among other factors, on temperature, 
moisture, amount of decomposer organisms and substrate quality, the 
latter including factors like percentages of lignin and nitrogen. The 
type of substrate, such as humus, dead organisms and various types of 
litters, also has a high influence. The position in the soil profile de- 
termines the amount of oxygen and carbon dioxide available, which both 
affect the decomposition rate. 


The temperature апа moisture influences are at present represented in the 
model by carbon functions (see sections 5.1.7 and 5.1.8), which will be 
modified in the light of field measurements. 


4.7.] Decomposition of litter 

Only root iitter is present in the FH-layer and both the origin and 

the diameter of the root litter affect its decomposition. Thus, the 
Calluna vulgaris are decomposed at about half the rate of Pinus sílves- 
tris roots (Berg 1977, in manuscript). Linaonberry roots and pine 

roots are for the purpose of the model assumed to have about the same 
decomposition rate.At present the overall decomposition rate is ap- 
proximated to 40% per year. 
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During the decomposition of a certain litter type the C/N quotients decrease 
from about 100 commonly found in the roots down to, perhaps, 50 as is found 
in the SOM (Staaf & Berg, in manuscript). This is due to the recycling of 
nitrogen in the sample, carried out by the microorganisms and the release 
of carbon as C0». The decomposition of a certain litter type thus starts 
initially at a high value of the C/N quotient, with nitrogen limiting. 


Under the conditions prevailing at Ivantjarnsheden, nitrogen is probably 

the limiting factor in all stages of decomposition. The C/N, quotient is 

used for the regulation mechanism. It has a higher numerical value (about 
90) than the C/N quotient. 


4.7.2 Decomposition of SOM or humus and dead organisms 

In the SOM the presence of proteins and certain parts of microbial litter 
give rise to a pool of easily degradable C-N compounds which are well mixed 
with less degradable substances like lignin and humus. The C/N quotient of 
the SOM as a whole is relatively constant and has been measured as varying 
around 50 (Staaf & Berg, in manuscript). 


When we regard the SOM as consisting of the separate substances, humus and 
dead organisms, the conditions for the compartments become different. Dead 
organisms have Ү(1)- and Y(2)-carbon as well as Y(3)-nitrogen and a C/N] 
quotient of 10-12. Humus contains almost only Y(2)-carbon and it is assumed 
that the amounts of protein-carbon and -nitrogen available are low. In the 
model they are given a low value which might be realistic and the resulting 
C/N4 quotient does not allow of any mineralisation. 


4.7.3 External breakdown 

In the breakdown of litter there will always be particles left which can 

no longer be called litter but are termed humus (see section 4.5). Me do not 
know at what level of decomposition a litter particle - a needle, say - 

will turn into nonidentifiable humus particles. The process is expressed 

as a fraction of the decomposition and is at present assumed to be 107. 


4.8 Carbon mineralisation (respiration) 

The carbon mineralised in the model is regarded as a fraction of the amount 
of carbon being decomposed. The fraction corresponds to the complement of 
the efficiency. If we write the efficiency as: 
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CASE SERE), the fraction thus is 1 - Bet for groth 

There are many values suggested for the efficiency factors for micro- 
organisms and values between 5% and 30% are given in the literature. 
According to Söderström (pers. comm.) a value of 10% could be expected 
for fungi, and has been used here mainly because the fungal biomass 
predominates (Bringmark et al., 1975). 


4.9 Decomposer biomass 

The biomass of the FH-layer at the test site is made up of microorganisms 
and soil animals, the latter with Acari and Cryptostigmata dominating in 
biomass, each varying between 300 and 20 ug/g soil. The microbial biomass, 
- fungi and bacteria - is dominated by the fungi (Bringmark e£ aZ. 1975). 
The fungal biomass as calculated by the Jones-Mollison technique (Jones 

& Mollison, 1948) varies between 6 and 24 mg: g7! of FH-layer and with 

1300 g of SOM we have 7-28 g:m"?. The living fungi appear to make up 
about only 10% of this amount (Söderström, pers. comm.). 


The biomass so far considered in the model is only the quantitatively 
dominating fungal one. The biomass variation is regulated by the growth 
and mortality functions. The whole biomass is considered to have a de- 
composing function in the system. 


4.10 Formation of decomposer biomass 

The growth rate of the microbial biomass is a function of the concentration 
of carbon and protein nitrogen and the microbial biomass. An any selected 
amount of microbial biomass and excess of carbon and nitrogen, the popula- 
tion will grow at a maximum rate. There is a dependence on temperature 

and moisture. At present these relationships are simple. See section 

5.3.7 and 5.1.8. 


4.10.1 Kinetics for substrate dependence 
At present two ways of expressing mathematically the kinetics of growth 
are used (see 5.1.10): 


- the Michaelis-Menten equation 
- the function for mean growth 


Both of these functions have been used by Parnas (1975a, 1975b) and the 
latter is a simplified version witnout the Michaelis-Menten parameters 
for carbon and nitrogen utilisation. Both approches modi- 

fy the value for maximum decomposer growth as environmentally adjusted 
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(5.1.9). They further take into consideration biomass and amounts of 
carbon and nitrogen available for decomposition. 


4.11 Mortality 


Proper description of the mortality process in a model is not possible 
at present owing to lack of knowledge and the following simplification 
is thus made: The microbial population has a constant average death rate 
under normal conditions (when an energy source is available). Thus, in 
normal conditions, death and growth occur simultaneously, and the net 
change in biomass depends on the relative magnitude of these two pro- 
cesses. When the populations are starved of energy - that is, carbon - 
the mortality constant increases and at the same time no growth occurs. 


4.12 Total nitrogen in soil water 

The total amount of nitrogen dissolved in soil water equals the amount 

of inorganic nitrogen dissolved viz. NH,* + N04, since the amount of 
organic N-compounds in solution is negligible (Bringmark & Petersson, 1975). 


The soil-water nitrogen varies considerably in amount from time to time 
depending,among other factors, on moisture conditions. 


In this model the soil water nitrogen is represented by the function Gl, 
which in model № will. Бе provided by module INORG. A very brief 
discussion follows here since details can be found in the forthcoming 
NIING report. 


fhe ionic concentrations can be determined for the soil system by assum- 
ing chemical equilibrium between the soil solution and the adsorbed phase. 
In this way, the amount of every ion (e.g., NH4") in soil water, will 

depend on the instantaneous state of the soil complex. 


There is a simpler, but only approximate,method of doing this. The varia- 
tions in soil-water concentration for a given ion can be determined start- 
ing from variations in the total soil concentration (e.g. as produced by 
mineralisation) by means of the soil buffer capacity. The soil buffer cap- 
acity for a given ion can in turn be determined from the adsorption iso- 
therm for the ion. 
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Work is currently in progress to obtain this kind of information concerning 
NH, and other ions. NIORG simply assumes that all їн in the model is in 
the soil-water phase (i.e. buffer capacity equal to unity and adsorbed 


amount equal to zero). 


4.13 Nitrogen mineralisation 

The transformation of organic nitrogen to ammonia takes place at a rela- 
tively low rate in the FH-layer at the Ivantjarnsheden site. The nitrifi- 
cation process has a rate much lower than that for formation of ammonia 
and denitrification does not take place at all. 


4.14 Transport of dissolved substances 

The amounts of carbon compounds and organic nitrogen transported between 
layers are negligible and are assumed to be equal to zero in the model. 
The NO.” and Анд? transport between the different layers is also assumed 
to be equal to zero in this model. Leaching of inorganic nitrogen will 
be developed in the NIING model. 


4.15 Immobilisation of inorganic nitrogen 

Immobilisation describes the rate at which inorganic nitrogen is 
transferred to microbial biomass. The immobilisation takes place only as 

a transfer from the soil-nitrogen compartment to the biomass compartment. 
It is assumed that the microorganisms can use organic and inorganic nitro- 
gen equally well but at rates depending on their relative availability. 


5. Mathematical description 


According to one of the hypotheses of this model, organic materials with 

a higher C/N, ratio than the critical value do not contribute to nitro- 
gen mineralisation. That means that, as far as nitrogen mineralisation 

is concerned, it is of no significance to have those materials as separate 
entities in the model or clustered together. Since the values of the 
parameters determining the critical C/N ratio are very uncertain, we 

have aimed at building a model with a rather flexible structure. Flexi- 
bility is a necessary condition if one wishes to test the significance 

of alternative structures in a model. 
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At its maximum capacity, NIORG can manage, at the chosen soil level, five 
organic materials (substrates) and five boxes for each material. The cho- 
sen structure for a given run can be specified by giving convenient values 
to certain parameters (see Table 1). If the number of boxes for a given 
material is set equal to zero, NIORG will cluster this material with SOM 
and the inputs corresponding to this material will correspondingly be 
attached to SOM. 


‘In order to have the desired flexibility in structure, it has proved neces- 
sary to attach only one transfer to each box. Thus, the model connects box 
Xnn to the transfer TOOnn which in that way actually becomes the time de- 
rivative of box Xnn. Each transfer (derivative) is defined as the alge- 
braic sum of a series of processes which are represented by the auxiliary 
functions G. We can understand this if we look at Fig. 1. The time deriva- 
tive of compartment NH,*, for example is calculated as 


T0001 = Z (619(m)-G20(m))-G70 where the sum extends over all substrates m, 


and in a similar way for the remaining compartments. This procedure is also 
convenient from the point of flexibility: processes can be changed by just 
redefining the corresponding G, or new processes can be included by just 
increasing the number of G-functions in the model. This has been the case 
with e.g. nitrogen fixation, nitrogen uptake by roots and transport between 
horizons which were not included at the preliminary version of the model. 


Functions Gl to G49 are reserved for the organic modules ORG and DECERS, 
G50 to G69 for INORG and G70 and onward for the ROOT module. A list of G- 
functions for NIORG is presented. The model needs some G-functions which 
should actually be calculated by some other module, e.g. Gl and G70. It 
should then be understood that these functions, as they are presented in 
NIORG, are very provisory and simplified im character. 


The modules (subroutines) around which NIORG is built are: 


- DECERS This subroutine calculates decomposers growth (biomass formation), 
mortality, etc. 

- ORG Calculates decomposition, mineralisation, inputs to the substrates 
organic transports, etc. 


The list of equations defining each G will mow be written with the para- 
meters corresponding to LITTER. ЈЕ is clear that each time it is called, 
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ORG fits the parameters according to the material under consideration. 

The same applies to compartments, i.e. they are fitted at each call of 
ORG; internally, ORG uses the five component vectors Y, with the following 
meaning: Y (1) protein-carbon for the given material, Y (2) other carbon, 
Y (3) protein nitrogen, Y (4) and Y (5) other possible chemical components 
for the given material. 


The model produces a rather large amount of information and it is probably 

| unnecessary to make a graph of each G-function. Even if all of them аге 
actually available for plotting, we have made a selection of what we con- 
sidered the most interesting information (see Lecend to NIORG graphs, p.46). 
The initial compartment values for a typical run are seen in Table 2 (only 
one soil level; in this case the FH-horizon). 


5.1 Auxiliary variables 


5.1.1 Inorganic nitrogen (NH,*) in soil water (g'm? soil) 
The inorganic nitrogen in soil water is here represented by the function. 


Gl. As mentioned previously, we simply have: 


G 1 = X01 
where: 
X01 is the total amount of inorganic nitrogen 


5.1.2 Total nitrogen available to decomposers (gm2 soil) 
G2L- Y3L + Gl 


where: 
Y3L is the protein nitrogen in litter 


5.1.3 Total carbon available to decomposers (gem? soil) 
G3L= YIL + Y2L 
where: 

YIL is the protein carbon in litter 

Y2L is other carbon in litter 


5.1.4 C/N, ratio in substrate 
CAL.= G3L/Y3L 


2] 


5.1.5 Total carbon used by decomposers per unit biomass increment 
G5 - P018/P019 
where: 
P018 is average carbon fraction in the decomposer 
P019 is ratio of carbon assimilated to carbon decomposed (assi- 


milation efficiency) 


5.1.6 Critical carbon-nitrogen ratio 
G6 = G5/P020 


where: 
P020 is average nitrogen fraction in the decomposer cell 


5.1.7 Temperature scale factor for decomposer growth 
It is given by function G7; to define this function, four parameters are 
required as shown in the figure and their values are given in table 1 


(see P005, P006, P007 and P008). 





-10 0 10 20 30 40 Temp. (°C) 
5.1.8 Moisture scale factor for decomposer growth 
G8_is given by the same kind of function as G7; the parameters in this 


case are: P009, P010. P011, P012. 


5.1.9 Maximal decomposer growth rate environmentally adjusted (аут 1) 
G9L- Р027 * G7 * G8 
Mhere: 

P027 is maximal (optimal) growth in LITTER 
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5.1.10 Decomposer growth (9 biomass:m ?.-day-l) 

There are at present two alternative ways of calculation included in N1ORG, 
and both of them are tested in different runs. If parameter P017 is set 
equal to zero, NIORG will calculate growth according to Parnas, (975a): 


- Decomposer growth (stationary) (g biomass · m? . day~!) 
GI3L = 691 * G3L * 621 * X02 / [(Р080+631) * (P085*G2L)] 
where: 


PO80 is the Michaelis-Menten constant for litter carbon utilisation 
P085 " "Ou Lo " n е nitrogen " 
X02 is decomposer biomass (g biomass -m?) 

If parameter P017 is set equal to unity NIORG will calculate growth 

according to Parnas (1975b): 

- Decomposer growth (mean) g biomass-m72-day~!) 

If G4L is greater than G6 the mean growth is calculated as: 

G13L = (G2L+G3L) / [(G4L+1) * G22L * p020] 

If G4 is less than G6 this growth is defined as: 

6131 = (G2L+G3L) * GAL /[(G4L+1) * G22L * G5] 

In these equation G22L is the growth time in litter (time taken for 

decomposers to reach the stationary phase) and is defined as: 


(621 +631) 


] Z G9L (days) 
5.1.11 Decomposer mortality (g biomass "72 :йаут1) 
If there is carbon available in the system: 


G12 = P032* X02 
where: 
P032 is normal mortality rate 


If there is no carbon available: 
G12 = P033 * X02 
where: 
P033 is starvation mortality rate 
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Mortality must go as input to DO, if this substrate is present in the 
chosen structure for the run, otherwise to SOM. To do that, 612 is con- 
veniently partitioned by the model into the different chemical components 
Y (1), Y(2) and y(3). 


5.1.12 Litter carbon decomposition (9-т-2- ау” 1) 
G14L = G5 * 6131 


5.1.13 Nitrogen fixation (g-m?.day^! 


G23L= P021 * P020 * G13L 
where: 


) 


P021 is the fraction of biomass X02 which fixes nitrogen 


5.1.14 Protein-carbon decomposition (g-m?.day^! 


) 
If G4Lis greater than G6 we have: 

6151 = (P034 * P020 * G13.* Y3L/G2L * (1 - P021) 
where: 


P034 is C/N ratio in proteins 


If G4L is less than G6 we have: 
G15L= 614% YIL/G3L 


5.1.15 Other carbon decomposition (g:m2-day*!) 
Gl6L= Gl4L- 6151. 


5.1.16 Respiration rate (g-m 2-day” 1) 
GI7L= (1 - P019) * Gl4L 


5.1.17 Protein-nitrogen decompostion (9-72. дау-1) 
G18L- 6151/ P034 


5.1.18 Protein-nítrogen mineralisation (g*m 2. day-l) 
6191= 6181- (P020 * G13L* Y3L/G2L) * (1 - P021) 


5.1.19 Nitrogen immobilisation (g*m'?-day"!) 
G20L- (P020 * 6131* G1/G2L)* (1 - po21? 
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5.1.20 Protein nitrogen net mineralisation (g*m?- даут 1) 
G21L- 6191- 6201 $ 


5.1.21 Total material (g:m ^) 
G30L- YIL + Y2L + ҮЗІ + Y4L + Y5L (at present Y4L and үбү are zero) 


5.1.22 Carbon-nitrogen ratio (measured) of substrate 
G32L= (G3L+ P018 * X02)/(Y3L + P020 * X02) 


5.1.23 External breakdown (Fragmentation) (g-m ^.day') 
We have assumed that this process can simply be expressed as a fraction 
of total carbon decomposition: 
G39 L= P037 * 6141 
where: 
P037 expresses this fraction 


Each component N of the substrate will accordingly be broken down as: 
G39L* YNL/G30L 


5.1.24 Transport of substrate between horizons (gm? day^!) 
Protein-nitrogen transport from S to FH horizon on a yearly basis is 
given by parameter P092. 


It has been assumed that this transport occurs only in the form of DO 

and at a constant rate during a period of approximately 200 days per year. 
Transport of other chemical components is calculated by knowing the C/N] 
ratio of the transported material (D0). 


5.1.25 Nitrogen uptake by roots (g:m?-day*!) 

By making certain assumptions, the ion concentration in solution at the 
root surface can be deduced from the ion concentration in solution in the 
soil surrounding the root (Baldwin et al., 1973). The concentration at 
the root surface is of primary interest if one assumes the rate of uptake 
as proportional to this concentration. 


In NIING, the ion concentrations in the soil solution are calculated, and 
this should be the basis for the development of the ROOT module. This 
module will also require a certain amount of parameters which are very 
uncertain at present. 
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For that reason, NIORG formulates nitrogen uptake by rocts in a very 
simple way; the main interest in this model is to observe the effect of 
this process on other processes in the model. 


We simply have: 


3 
G70 = P077 * X б19 (m) 
m=1 


where: 
P077 represents the fraction of mineralisation which goes to roots 
and 619 (т) is the mineralisation from substrate m. 


By varying P077 between zero and unity. we can observe the effect of uptake 
in decomposer growth, mineralisation etc. 


5.2 Driving functions 

As can be observed from the previous paragraph, in order to define G7 and 
G8 we need temperature and moisture at the given soil horizon. These, 
together with LITTER input are the three driving functions for NIORG. 
They should be provided by other models, but for the time being we have 
simulated them with very simple analytical functions of time. These are 
illustrated in Fig. 2: LITTER input (formation) at FH-layer in glitter’ 
m7 2+day~! (Fig. 2a), TEMPERATURE at FH-layer in © С (Fig. 2b) and 
MOISTURE at FH-layer in 9-97! (Fig. 2c). The time axis in these figures 
covers a period of two years. Several parameters are needed in order to 
define these functions. Some of them can be found in ТаЫе 1: P059, P060 
P065, P066, P071 and P072. Further, some parameters are needed to enable 
the partitioning of LITTER input into different chemical components. 
They can also be found in Tablel: P050 and P051. 


6. Model performance 


The model has been run a number of times and we have selected five of 
the runs for presentationand discussion here. 
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The first run (number 54) is the standard run, i.e. the гип carried 
out with the values of parameters and compartments presented in tables 
1 and 2 respectively. The sensitivity of the model to changes in 
different parameters and to changes in structure has been the object of 
the runs 55-58. We have investigated the sensitivity of the model 
(especially of mineralisation) to: 
- changing the way of calculating decomposers growth (run number 55) 
- сһапд1пд the amount of nitrogen uptake by the roots (run number 56) 
- changing the decomposers' assimilation efficiency (run number 57) 
= changing the structure of the model by combining DO and HUMUS in- 

to SOM (run number 58). 


At runs 54 to 57 three substrates - LITTER, DO and HUMUS - each with 
three compartments are present; that means P003 = P004 = 3. 


The outputs from each run of the model are shown in figures 3 to 7. 

The meaning of each variable symbol appearing in the figures is 

explained in legends to the graphs. Each run covers a period of two years, 
starting at July 1 ' of a year. Euler integration has been used and the 
time step has been set equal to one day. The execution time has been of 
the order of three minutes for each run. 


In order to clarify and organizé the discussions around the numerical 
information provided by the model, two numerical flow diagrams are 
provided for each run. One flow diagram is for carbon, the other for 
nitrogen. The figures in each diagram give the accumulated amount of a 
given process during the two year. period of the run and the unit is 
gm? per two years. 
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The accumulated processes for the carbon diagram are: 










carbon 
respired 







Total carbon 
in biomass 
formed 


For the nitrogen diagram we have: 
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6.1 Standard гип - run number 54 
In this run we have taken: 


- P017 = 1., i.e., biomass formation is calculated with the function for 
mean growth (5.1.10) and 
- P077 = 0., i.e., nitrogen uptake by roots has been set equal to zero. 


We can see from fig. 3-a that the model appears relatively stable, at any 
rate during the two year ‚гип. The steady increase in the compartment for 
inorganic nitrogen (X01) is caused by the absence of uptake by roots; 

670 being zero all the time, Fig. 3-b. 


The behaviour of compartment (X02) is not in total agreement with known 
facts from the Jaádraás site. Fungal biomass, which predominates, is low in 
the dry period in July and August and increases in the autumn and in April- 
May, that is, the biomass is high during the wet periods (Söderström, pers. 
comm. ). 


It is difficult to comment on the compartments for litter as we do not yet 

have information about the root litter dynamics. The compartments for dead 

organisms and humus show little variation and it is difficult to verify the 
variation obtained. 


The remaining G-functions have a typical form (Fig. 3-b) caused by the tem- 
perature and moisture modulation of biomass growth. 


In Fig. 3-c we can observe several facts which are direct consequences of 
the hypoteheses in the model. 


With the values given to parameters P018, P019 and P020 (Table 1), the re- 
sulting critical C/N ratio is equal to 100. The initial C/N] values for 
LITTER, DO and HUMUS are 90, 10 and 3204 respectively. This means that only 
LITTER and DO will mineralise nitrogen; which is confirmed by the plots 

of GI9L, 61900 and G19H. One of the conclusions presented by Parnas (1975a)* 
is that a material with a lower carbon-nitrogen ratio than the critical 
value will not change its carbon-nitrogen ratio as decomposition 


* Most of the conclusions presented here are of a general nature, regard- 
less of the actual way used to calculate decomposers growth. They are 
a consequence of tne hypothesis on which the model is based. So, such 
conclusions are valid here, even if we at this run have used mean growth 
instead of stationary growth to calculate biomass formation. 
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goes on. This is confirmed by the plots of G4L and б4р0, On the other 
hand, a material with a higher carbon-nitrogen ratio than the critical 
value will have an increase in its carbon-nitrogen ratio if none or 
little NH,” is present in the system. This is confirmed by the plot of 
G4H (HUMUS carbon-nitrogen ratio): it increased from 3204 (initial) to 
3208 (final). This has been so, even considering the break-down of LIT- 
TER and D.0. which should tend to decrease the HUMUS carbon-nitrogen 
ratio. 


The "measured" C/N ratio (substrate plus biomass) behaves in a more 
complicated way: 6321, 632900, G32H The variation of the C/N ratio in 
litter (G32L) appears reasonable. 

The decomposer time of orowth (5.1.10. ) was more than two months 
for the three substrates LITTER, 0.0. and SOM-carbons (G22L, 62200 
G22H) - 


The carbon and nitrogen diagrams for this run are as follows: 


90.74 


81.64 


9.10 





The contributions from LITTER,D.0. and HUMUS to the total decomposed 
carbon (90.74) were 76.22, 14.34 and 0.18 respectively. Respir- 
ation and mineralisation mean values for the period of the run are 
0.1134 and 0.0023 g - m? day ! respectively. The respiration value 
(83 g for 2 years) might be low compared with the value of 120-200 g C 
suggested by Staaf (pers. comm.). Unfortunately it is not possible 
to evaluate the mineralisation rate. 
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Conclusion 4 by Parnas (1975 a) states: "The rate of mineralisation is 
increased as C/N, is decreased". This is confirmed here: To the 

1.66 g m of nitrogen mineralised during the two years the contri- 
bution from LITTER and D.0. were 0.33 and 1.33 g m"? respectively. 


6.2 "Stationary growth" for decomposers - run number 55 
In this run we have 


- P017 = 0. i.e. biomass formation is calculated with stationary 
growth (5.1.10. ), and 
- P077 = 0. as before 


Fig. 4-a shows that the model behaviour is not as good as in the pre- 
vious run: LITTER and D.0. compartments increase with time, while 
BIOMASS compartment (X02) decreases. 


With the exception of that, most of what has been said about general 
trends for the previous run also holds for this one. 


The carbon and nitrogen diagrams are now: 





The contributions of each material to carbon are: 20.54 for litter, 
7.38 for D.0. and 0.31 for humus. Regarding nitrogen mineralisation 
the contributions are: 0.68 for D.0. and 0.06 for litter. The respir- 
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ation and mineralisation mean values are 0.040 and 0.001 g m ^ day 


respectively. 
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The lower figures from this run are a consequence of the behaviour 
of the compartments. One possible way to improve this could be to 
change the Michaelis-Menten parameters (see Table 1), the mor- 
tality parameter or the optimal growth parameters. Besides the com- 
plications involved, it is known that the use of this kind of formu- 
lation to simulate decomposer growth in natural systems has been 
criticised.Thus, at this point, a sound decision seems to go on with 
the formulation of growth used in run 54. In the runs to follow it 
will be understood that the mean growth is used to calculate 
biomass formation. 


6.3 Effect of nutrient uptake by roots - run number 56 
Now, we have: 


- P077 = 0.9 (uptake by roots is included) 


The effects of this change in the root parameter can be seen in 
Fig. 5-a; the NH4 
er than in run 54. The nitrogen uptake rate is shown by G70 in 
Fig. 5-b. 


* values (X01) are now one order of magnitude low- 


The general trends discussed in run 54 are, of course, not modified 
by the new process and decreased NH,* - level. 


Some new, interesting facts can be noticed, however. Before these 
are discussed, we present the flow diagrams: 


90.55 


81.45 ' 


9.10 
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As before, the contributions to decomposed carbon are: 76.08 for litter, 
14.30 for 0.0. and 0.17 for humus. The contributions to mineralisation 
amount to 0.16 for litter and 1.30 for р.0. Respiration and mineralis- 
ation mean values are 0.1132 and 0.0020 g me day ! and thus, respiration 
mightalso be low here with 83 g versus 120-200 as estimated by Staaf 
(pers. comm.). 


One major change that has taken place in the nitrogen flow, is that im- 
mobilisation is considerably lower than it was in run 54: 0.07 against 
0.27. This is not surprising considering that NH,* has been removed 
from the system. Since the rate of immobilisation is set to be pro- 
portional to the instantaneous fraction of NH,” in the system (5.1.19) 
the result is not surprising. Less expected is the decrease in mineral- 
isation rate: 1.46 against 1.66 at run 54. This is, however, in agree- 
ment with one of the conclusions given by Parnas (1975 a). This is: 
"Conclusion 5. The higher the amount of NH,” in the area, the higher 
will be the rate of mineralisation". In other words: the higher the 
amount of nitrogen taken up by the root system, the lower will be the 
rate of nitrogen mineralisation. 


It is important to notice that the biomass dynamics has practically not 
been influenced by the new process: the mean value of X02 has been 
16.95 for run 54 and 16.94 for the present run. The biomass formed has 
been 18.15 at run 54 and 18.11 at the present run. 


This has two important consequences: 


- first, the amount of nitrogen required for growth is the same at both 
runs: 0.91. As immobilisation is now lower, the decomposers must have 
recourse to organic nitrogen: 0.82 against 0.62 at run 54, 

- second, the decomposed amount of organic nitrogen is the same for 
both runs: 2.28. As a consequence of all this; the organic nitrogen 
surplus turns out to be lower at this run. 


6.4 Effect of changing decomposer assimilation efficiency - run 

number 57 
One important parameter in the model is the biomass efficiency in car- 
bon utilisation: P019. It derives its importance from, amonq other things, 
the fact that it is one of the parameters determining the critical C/N 
ratio (5.1.6). 
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The higher this parameter value, the lower will be the critical C/N 
ratio in the model. This could indicate that the model should be 
very sensitive to changes in this parameter. 


In order to test this, P019 was doubled in this run and thus the 
critical carbon-nitrogen ratio went down to 50.Thus,in this run: 


- P019 = 0.2 

Fig. 6-c shows some expected results: 

- 6191 equals zero; i.e. there is no mineralisation from litter. 
This lack of mineralisation from litter might differ from reality 
(Popovic, 1976) 


- GAL (C/N, in litter) increases from 90 (initial) to 96 


Fig. 6-a shows that the litter compartment tends to increase at 
this run. 


The flow diagrams are now: 


36.41 1.25 


2913 





7.28 


E 






The contributions to total carbon decomposed are: 28.46 from litter, 
7.88 from 0.0. and 0.07 from humus. Clearly, in this case all miner- 
alisation came from 0.0. The mean values for respiration and miner- 

alisation are 0.040 and 0.001 g ar day”! respectively. The respir- 

ation value (30 g for 2 years) apparently is low compared with that 

given by Staaf (pers. comm.) on 120-200. 
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6.5 Effect of simplifying model structure; HUMUS and D.0. combined in 
SOM - run number 58 

As stated previously, in this run 0.0. and humus are clus- 

tered together in SOM. 


Fig. 7-a shows that the litter compartment behaved, approximately in 
the same way as in run 54. The SOM compartment increased with time; 
this is due to the fact that the value for optimal growth on humus 
has now been used for SOM (a mean value between those of D.0. and hu- 
mus should, instead, have been used). 


The fact that mortality now actually gave an input to SOM by the model can be 
verified by looking at graph 845 in Fig. 7-с. It is due to this in- 

put that the carbon-nitrogen ratio in SOM decreases from 466 (in- 

itial) to 356 (final). 


Mineralisation from SOM - 6195 - is equal to zero, as it should be. 
The flow diagrams for this run are: 


75.45 


67.01 


7.44 





To the total carbon decomposed, LITTER has contributed with 

73.11 g m ^ and SOM with 1.34 g m^. All mineralised nitrogen now 
came from litter. This appa.ently is not true for the FH-layer 
(Popovič, 1976). The respiration and mineralisation yearly mean 
values were: 0.0930 and 0.0002 g її © day” | respectively - the respir- 


ation might be low. 
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7. Conclusions 
The runs discussed in the previous section show that: 


- the model behaves, at present, in the way that should be expected 
from theoretical (analytical) considerations 

- the model is a convenient tool to analyse the effects produced on 
some processes when changes in parameters and in structure are 
introduced 

- the model predicts values for decomposers respiration and nitrogen 
mineralisation which seem to be lower than estimated values obtained 
from measurements. 


At this point, what seems to be needed most are experimental data to 
further modify, build, and verify the model. For planning this further 
work a convenient next step would be to carry out a more extensive 
sensitivity analysis with the present model. 
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Comments on some parameters used (Table 1) 


There are too few measurements on processes and compartments from the 
Ih VA site to allowof too deep an analysis of the values obtained in 
the runs. A short discussion of some parameter values might be 
fruitful. | 


One important parameter is the assimilation efficiency (P019). This 
parameter is unfortunately very uncertain. The fractions of carbon 
and nitrogen in decomposer biomass (P018 and P020) are difficult to 
determine and it is possible that the literature values obtained 
might have to do. The optimum growth values in litter, humus and on 
dead organisms (P027-P029) are not determined - the values are chosen 
by making a rough estimate. The mortality rates (P032 and P033) are 
taken from Parnas and Redford (1974). The work by Söderström 

on estimating the death rate for fungi might give or confirm the 
size of this process. 


P037 is estimated from some preliminary decomposition values and P038 
is just a rough estimate. | | sx 


PO50-P051 in root litter are nct determined and needle litter values 
are used here. 


The relative response of microorganisms to moisture and temperature 
(PO005-PO12) are so far represented by two simple functions. This re- 
sponse has a central role in the model. 


The dynamics of the root litter input is shown in Fig. 2 a. The input 
values (P059, P060) are given asa first approximate value by Persson (pers. comm. 
and the input dynamics are so far approximated due to lack of infor- 
mation. This parameter will be improved shortly. 


The fraction of the microorganisms that can fix nitrogen (P021) is 
setat 2%, a figure thatwill probably do. 


The transport of solid nitrogen from the S- to the FH-layer (P092) is 
discussed in section 5.1.24. The calculation is made from incomplete 
data but the value will hopefully be improved shortly. 
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Table 1. The most important parameters in Ї{10Вб 


Symbol 


P001 
P002 


P003 
P004 
P005 
P006 
P007 
P008 
P009 
P010 
P011 
P012 
P017 
P018 


P019 
P020 


P021 


P027 
P028 
P029 
P032 
P033 
P034 
P037 
P038 
P050 


P051 


P059 


Meaning 


No. of litter fractions 


Code parameter to define 
structure 


No. of litter or humus boxes 


No. of D.0. boxes 


Temperature scale factor (5.1.7) 


Switch (see 5.1.10) 


Fraction carbon in biomass 
(5.1.5) 


Assimilation efficiency (5.1.5) 
Fraction nitrogen in biomass 


(5.1.6) 


Fraction biomass which fixes 


nitrogen (5.1.13) 


Optimal growth in litter(5.1.9) 


Be 0.0. 


humus 


Normal mortality rate (5.1.11) 


Starvation rate (5.1.11) 
Protein C/N ratio (5.1.14) 
Litter breakdown (5.1.23) 
0.0. п п 


Fraction protein carbon in 
litter input at FH level 


Fraction other carbon in 

litter input at FH level 

Daily litter input during 
spring (FH level) 


Typical References 
value 


1 


„9 Söderström, pers. comm. 


0.1 "s 
0.05 s 


0.02 Granhall, pers. comm. 


„005 day”! 

.0015 даут! 

.0001 day"! 

.0015 day^lParnas & Redford (1974) 
.02 days! — "- 

Parnas (1975a) 


O O O A oocoosZono 


1.5 g т^”2дау^! Persson, pers. comm. 


Symbol 


P060 
P065 


P066 
P071 
P072 
P077 


P080 


P081 
P082 
P085 


P086 
P087 
P092 
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Meaning 


Daily litter input during 
autumn (FH level) 
Maximum temperature of year 


Minimum temperature of year 
Maximum moisture of year 
Minimum moisture of year 


Mineralisation fraction which 
goes to roots (5.1.25) 


Litter carbon Michaelis-Menten 
parameter (5.1.10) 


0.0. "- 
Humus "- 


Litter nitrogen Michaelis-Menten 
parameter (5.1.10) 


D.O. "= 
Humus "- 


Protein-nitrogen transport from 
S to FH horizon (5.1.24) 


Typical References 
value 


1.5 g m72day7! Persson, pers. comm. 


15.9C Measured value 
(Problem area 2) 
-3.°C "- 
4 ъ -1 "o 
- g'g 
0.5 g:g^! "- 
D to 1 


20.0 g'm? Parnas (1975a) 


1.0 g-m? = 
20. °т° "o 
1.0 g:m-? "- 


1.0 gm"? "- 
3 -2 "o 
1.0 gem 


0.145 g m. угт! present paper 
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Comments on the initial compartment values (Table 2) 


It appears that in the FH-layer we have measured initial values for the 
compartments fairly well - an exception is of course the compartment 
for dead organisms in which case the correct order of magnitude would 
be as close as we could hope to come. The dynamics of the compartments 
is less well known. 


Inorganic nitrogen in the FH-layer has been measured by Popovic (1975, 
1976) and varies between 27 and 53 mg/m* in the period from May and to 
November. : 


The litter compartments (YIL,Y2L,Y3L) have been measured ` (Persson, pers. 
comm.) and the values obtained appear to do for the model. The dynamics 
of the root litter compartment which will be determined will be very 
useful to verify the variation obtained by the model. There is a lack 
of analysis on root litter, however, and the values for protein-carbon 
and-nitrogen as well as other carbon were taken from analysis carried 
out on needle litter. It is possible that they are not correct. 


The biomass compartment is dominated by fungal and bacterial biomass 
(Bringmark e£ aZ. 1975) and only the values for these are used. The 
variations during the year might be large- it is known that the biomass 
decreases in the dry period from June to August and is high in spring 
and in autumn (Söderström, pers. comm.). Owing to lack of information the 
size of the compartment has, as an approximation, been given the initial 
value of double the fungal biomass. 


The humus compartment has been measured (Staaf & Berg, in manuscript) 
and it appears that the value obtained is good enough. 


The compartment for dead organisms has been roughly estimated. as a 
fraction of the biomass compartment. 
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Table 2 NIORG compartments (initial values) 
The compartments and their initial values for a typical run. The struc- 


ture for this run is: one soil horizon (FH), three organic materials, and three 
boxes for each material. 


Symbol Meaning Value Reference 

X01 NH] 0.1 gm Popovic (1975, 1976) 
x02 Decomposer biomass 15.0 *- : Söderström, pers. comm. 
YIL Litter protein-carbon 7.0 te present paper 

Y2L Litter other carbon. 150.5 "- "- 

Y3L Litter protein-nitrogen 1.75 "- "- 

Y1DO D.0. protein-carbon 12.0 "- Soderstrom, pers. comm. 
Y2D0 D.0. other carbon 18.0 “a "- 

Y3DO D.0. protein-nitrogen 3.0 > "- 

YIH Humus protein-carbon 2,0 "a present paper 

Y2H Humus other carbon 1600.0 "- "- 


Y3H Humus protein nitrogen 0.5  "- 5s 
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Figs. 3-7 Legends to NIORG graphs 


Symbol 


613Т 
617Т 
G19T 
G20T 
G21T 
G23T 
G70 
6131. 
6221. 
G14L 
G30L 
GAL 
6321. 
Gl 9L 
61300 
62200 
61400 
63000 
G4D0 
G32D0 
61900 
G13H 
G22H 
СТАН 
G30H 
G4H 
G32H 
G19H 
X01-Y3H 


Meaning 


Total biomass growth 

Total respiration 

Total nitrogen mineralisation 
Total nitrogen immobilisation 

Net nitrogen mineralisation 
Total nitrogen fixation 

Nitrogen uptake by roots 

Biomass growth from litter 

Growth time in litter 

Litter carbon decomposition 

Total amount of litter 

Litter C/N, ratio 

Litter C/N ratio (measured) 
Nitrogen mineralization from litter 
Biomass growth from 0.0. 

Growth time in 0. 0. 

0.0. carbon decomposition 

Total amount of D.0O. 

D. 0. C/N, ratio 

0.0. C/N ratio (measured) 

Nitrogen mineralisation from D.0. 
Biomass growth from humus 

Growth time in humus 

Humus carbon decomposition 

Total amount of humus 

Humus C/N, ratio 

Humus C/N ratio (measured) 
Nitrogen mineralization from humus 
Compartments as defined in table 2 


Units 


biomass m ^ day”! 


9 

g carbon m? day 

g nitrogen me day 
n 

g biomass m? day”! 

months 

g carbon m? day ! 

дет” 

g nitrogen т^? day 

g biomass m? day ! 

months 

g carbon m? day! 

д-т 

g nitrogen m"? дау” 

g biomass m^ day ! 

months 

g carbon m? дау”! 

g: m" 

g nitrogen m? дау” 


1 


1 


1 


1 
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Fig. 4-b Variables defined in legends:Figure to the right of the variable symbol gives 
the value of the variable at the top of the axis. Time starts at July ] ofa 
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Fig. 5-b Variables defined in legends. Figure to the right of the variable symbol gives 
the value of the variable at the top of the axis. Time starts at July 1: ofa 
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Fig. 5-с Variables defined in legend, Figure to the right of the variable symbol gives 
the value of the variable at the top of the axis. Time starts at July 1 -ofa 
year and ends two years later; each step оп the time axis represents 200 days. 
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Fig. 6-b Variables defined in legends.Figure to the right of the variable symbol gives 
the value of the variable at the top of the axis. Time starts at July 1 -of a 
year and ends two years later; each step . on the time axis represents 200 days. 
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Fig. 6-c Variables defined in legends.Figure to the right of the variable symbol gives 
the value of the variable at the top of the axis. Time starts at July 1 «ofa 
year and ends two years later; each ‹ Step on the time axis represents 200 days. 
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Fig. 7-b Variables defined in legends. Figure to the right of the variable symbol qives 
the value of the variable at the top of the axis. Time starts at July l- ' of a 
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